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OA ' Abstract. Stars between about 4 and 25 solar masses spend a significant 

$_i \ fraction of their post-main sequence lifetime as red supergiants (RSGs) and lose 

C^H' material via stellar winds during this period. For RSGs more massive than 10 

' solar masses, this mass loss becomes of evolutionary significance, and probably 

determines the upper mass limit of RSGs in the Hertzsprung-Russell diagram. 
Despite decades of observations, the driving mechanism responsible for mass loss 
in RSGs remains unknown. Mainly this is because the optical spectrum accessi- 
ble from the ground provides almost no useful wind diagnostics, and what infor- 
py^ mation is obtained is spatially averaged over the entire wind volume. However, 

within the last decade, Hubble Space Telescope (HST) observations of many use- 
ful ultraviolet wind diagnostics have been obtained at a high signal-to-noise ratio 
and spectral resolution. In particular, RSGs in eclipsing binaries can provide 
spatially resolved observations of stellar chromospheres and winds. I review pos- 
q . sible RSG wind acceleration mechanisms, discuss some observational constraints, 

' and present some empirical models of RSG chromospheres and winds. 

1. Stellar Winds in the H-R Diagram 

m ; 

\ Stellar winds are a ubiquitous phenomenon across much of the Hertzsprung- 

Russell diagram (HRD). These winds can be broadly grouped into the three 
main categories of hot star winds, coronal (solar-type) winds, and the cool winds 
Tj- ■ of evolved, late- type stars (Figure [T]). 

o . 

1.1. Hot Star Winds 

Stars earlier than spectral class B3 typically have high velocity winds, and sig- 
nificant mass-loss rates. Radiation pressure is substantial in these stars since 
the (large) opacity of iron-group metals peaks in the mid-ultraviolet where the 
radiation field is strong, and so Py»A is th e dominant term in the momentum 
equation. ICastor. Abbott. &: Klein (|l975l ) showed that the winds of hot stars 



are driven by the force of radiation on optically thick atomic lines in regions 
above the photosphere, where the continuum is optically thin. 

1.2. Coronal (Solar-type) Winds 

Main sequence stars later than about F0 have substantial convective zones, 
which interact with stellar rotation to produce magnetically active photospheres. 
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Figure 1. Stellar Winds in the Hertzsprung- Russell Diagram. 



Acoustic waves and turbulence generated by convective "noise" propagate out- 
ward into the chromosphere, where the acoustic waves steepen into shocks and 
dissipate, producing localized, transient heating of the chromospheric gas (and 
the non-magnetic chromospheric line emission). Gas in the corona, lying above 
the chromosphere, is heated to million degree temperatures by dissipation of 
magnetic flux and reconnection of magnetic flux tubes tied by their photospheric 
footprint to convective motions. Coronal winds are tenuous with a mass-loss rate 
too small to be of evolutionary importance. These are driy en mainly by the ther- 
mal pressure gradient of the hot corona, the Parker ( 19581 ) mechanism, although 
dynamic events in the corona (magnetic reconnection events, flares, and coronal 
mass ejections) also contribute to the solar wind mass flux (and presumably to 
that of similar stars also). 



1.3. Cool Winds of Late- Type Stars 

The red supergiant binary a Her was first shown, by iDeutschi ( 19561 ). to be 
losing mass via a stellar wind. The intervening decades of research have failed 
to produce a viable ph ysical process that exp lains the mass-loss process in red 
giants and supergiants ( Judge Sz Stencellll99ll ). In these stars, mass loss occurs 
via a low velocity (typically ~20-50 km s _1 ), massive (10~ 9 to 10~ 5 M Q yr _1 ) 
wind. This cool wind regime i s sep arated in the HRD from that of coronal 
winds by the iLinskv Sz Haisch (1979) dividing line (Figure [1]), which lies near 
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spectral class Kl for red giants and around FO for the most luminous supergiants. 
Noncoronal stars show generally weaker chromospheric emission line fluxes and a 
lower degree of ionization than the coronal stars, and have little or no detectable 
X-ray flux. There is little evidence of emission from gas warmer than about 
20,000 K, and most of the gas must be much cooler than this. 



2. Possible Mass-Loss Mechanisms in Red Giants and Supergiants 

The most stringent constraint on possible mass-loss mechanisms is provided by 
the slowly accelerating nature of the wind. Consequently, most of the work 
done in driving these cool star winds goes into lifting material out of the stellar 
potential and not into the kinetic energy of the ejecta. This implies an effective 
driving force that only slightly exceeds gravity over distances of several stellar 
radii. 



2.1. Pulsation 



Pulsation plays a dominant role in initiating mass loss in Mira variables, OH/IR 
stars and other cool, evolved stars with optically thick, dusty winds ([Bowen 



1988). These are mostly highly evolved asymptotic giant branch (AGB) stars, 



but probabl y include some massive red supergiant (RSG) stars also, such as 
WOH G64. IVan Loon et al.l ( 20081 ) showed that the ratio of wind kinetic en- 



ergy to the power supplied by pulsation (^Mv^/ '-E pu is) varies smoothly with 
pulsation period, being uniformly large for periods exceeding 600 d. The RSGs 
blend in indistinguishably from the AGBs in this trend, suggesting that pulsa- 
tion may be important in RSGs also. This is consistent with the suggestion of 
Judge Sz Stencell (|l99ll ) that long period waves (with periods longer than the 



acoustic cutoff period) g enerated by ph o tosph eric variability may drive mass 
loss. On the other hand. Ijosselin Sz Pled ( 20071 ) state "the extrapolation of the 



theory of mass loss of AGB stars [to RSGs] seems irrelevant as RSGs have only 
irregular, small-amplitude variations and are not pulsating in a similar manner" . 
In summary, while it is clear that pulsation is a leading contender amongst pos- 
sible wind-driving mechanisms, and it is certainly energetically viable, there still 
is no definitive observational evidence that the winds of (most) RSGs are actu- 
ally initiated or driven by pulsation. Nevertheless, the observational record is 
sparse, and the role of pulsation in RSGs remains very much an open question. 



2.2. Radiation Pressure on Dust Grains 

Dust grain opacities, unlike atomic opacities, peak in the infrared, and so en- 
sure good coupling of the radiation field to the circumstellar material. As such, 
dust-driven winds (usually initiated by pulsation) have been among the most 
popular of proposed mechanisms exp l aining mass lo ss in red giants and super- 
giants, e.g. Lafon Sz Berruver ( 199 ll ). iHofnerl ( 20081 ). There is much evidence 



that dust-driven winds are the operative mechanism in Miras, OH/IR stars, and 
the optically thick, dusty winds of other highly evolved AGB stars (and possibly 
a few RSGs such as WOH G64). But these conclusions can not possibly be 
generalized to the very different conditions found in the optically thin winds of 
RSGs (and first ascent red giant branch stars). Some RSGs have relatively little 
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dust, and where the dust emission has been sp atially resolved by IR interferom- 
etry in the case of the M2 supergiant a Ori (Be ster et al.|[l99a ). the observed 



dust shell lies at a dista nce of ~ 30.R*, far exterior to the empirically inferred 
wind acceleration region ( Harper. Brown. Sz LimlbOQll ). Furthermore, VV Cep, 



another M2 supergiant in a binary system and spectral proxy of a Ori, contains 
little or no dust, yet has a comparable mass-loss r a te to that of a Ori. Similar 
argum ents were advanced years ago by iGoldberd (|l979l ) and Uudge fc Stencell 



(1991), and the sassage of time has not diminished their logic. Radiation pres- 



sure on dust grains can not play a significant role in driving the optically thin 
winds of RSGs. 

2.3. Dissipation of Alfven Waves 

One of the most plausible mechanisms proposed to drive mass loss in RSGs 
has been the dissipation of Alfven waves. Unlike acoustic waves, Alfven waves 
can have large damping lengths that are far more effective in transferring en- 
ergy and momentum to the circumstellar gas over distances of the order of a 
stellar radius, as needed to pro duce stellar winds consistent with observations 
(jHartmann fc MacGregorlll980l ). Alfven wave dissipation is predicted to result in 



extended, warm chromospheres in l ow gravity stars — a result t hat was thought 
to be consistent with observation. Hartmann &: Avrettl ( 19841 ) constructed an 



Alfven wave model appropriate for a Ori which predicts a warm wind with elec- 
tron temperatures peaking at T e ~ 8000 K at AR* (i?* is the radius of the stellar 
photosphere), and remaining above 5000 K out to 10-R* . However, th e se res ults 
are in strong disagreement with the empirical model of Harper et al.l (|200ll ) , in 



which the wind reaches a maximum temperature of just 4000 K in the lower 
chromosphere, then falls to 2500 K at 3R* and to 800 K by 10-R*, with a nearly 
constant ionization fraction of n e /nu ~ 10 -35 . It seems unlikely that Alfven 
waves, which couple only to ionized gas, could transfer the necessary energy and 
momentum to the wind in a nearly neutral gas that must be highly damping. 

2.4. Acoustic Waves 

Photospheric convection in late-type stars generates acoustic turbul ence, and 



some of this acoustic energy propagates upward into the chromosphere. lAthav &: White 



(1979) showed, from solar UV observations, that the solar aco ustic wave flux is 



barel y) sufficient to heat the chromosphere, but not the corona. lHartmann fc MacGregor 



(1980) computed acoustic flux models for low gravity stars and concluded that 
"the primary dynamical effect of propagating acoustic waves and the shocks 
that develop from them is to extend the essentially h ydrosta t ic por tions of the 
lower atmosphere, rather than drive a steady wind". ICuntzl ( 19901 ) confirmed, 



using a 1-D hydrodynamic code, that short period acoustic waves dissipate and 
heat the chromosphere but do not contribute significantly to mass loss. More 
sophisticated hyd r odyam ical modeling of the solar chromosphere carried out by 
Carlsson Sz Steinl ( 19951 ) further confirmed this outcome. The conclusion: short 



period acoustic waves heat chromospheres but do not cause significant mass loss. 

2.5. Radiation Pressure on Atomic and Molecular Lines 

The coupling of the cool radiation field of late-type stars to the line opacity of 
circumstellar gas is poor, and so the resulting radiative force is small. Molecules, 
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with strong infrared bands, such as carbon monoxide, provide a better match, 
but their overall abundance is small, and hence their contribution to the radiative 
force is limited. Typically, empirical models of the outer atmospheres of RSGs 
give a radiative acceleration of the stellar gas of a few times 10~ 2 cm s -2 , or 
loga ra d ~ —2. This is 1.5 dex below the log g of the lowest gravity supergiant. 
Radiation pressure on the gas alone can not drive a stellar wind in red giants or 
supergiants. 



3. A Tale of Four Stars 

Ultimately, the nature of the driving force accelerating cool stellar winds will be 
resolved by observation. To this end, I have analyzed ultraviolet observations of 
three bright red supergiants observed with the HST Goddard High Resolution 
Spectrograph (GHRS) and the Space Telescope Imaging Spectrograp h (STIS). 
With these results, and including the analysis of lHarper et al.l ( 200ll ) for a Ori, 



we can present comparative empirical results for the following four RSGs: 

• A Vel - the brightest K supergiant in the sky: 

V = 2.21 mag, K4 lb 

• a Ori - the brightest M supergiant in the sky: 

V = 0.50 mag, M2 lab 

• C Aur - the brightest eclipsing K supergiant binary: 

V = 3.75 mag, K4 lb + B5 V, P = 2.66 yr 

• VV Cep - the brightest eclipsing M supergiant binary: 

V = 4.91 mag, M2 lab + B, P = 20.34 yr 

The binaries are of particular interest because, in principle, they provide the 
opportunity of obtaining spatially resolved observations of the outer atmospheres 
and winds of these red supergiants. 

3.1. A Velorum 

The UV spectrum of this red supergiant, observed in 1994 by HST/GHRS, 
contains numerous emission lines of singly ionized metals, especially Fen. Op- 
tically thin lines have a pure emission profile, but optically thicker lines show 
self-reversed absorption. This absorption feature forms just shortward of the 
line centre for slightly optically thick lines, and deepens and shifts blueward 
with increasing line optical depth. We interpret the velocity at the blueward 
edge of the absorption trough to be the velocity of the wind where the line be- 
comes optically thin and the underlying emission profile (which forms close to 
the star) starts to show through. From this, we directly infer v(Nn), where v is 
the wind velocity and iVn is the hydrogen column density of the wind. Then, 
assuming a spherical steady wind with radial hydrogen column density Nn(r) 
and velocity v(r): 

T ~ 4ir»m H tf* {r) v(N H )dN H 1 

and so r(v) is found. Here \i is the mean particle mass, and tbh is the hydrogen 
atom mass. Inverting this result yields the desired empirical wind velocity law 
v(r) with a mass-loss rate M = 1.1 x 10 -9 M Q yr -1 and a terminal velocity of 
about 30 km s _1 (Bennett 2010, in preparation). 
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3.2. a Orionis 



The line profile analysis technique used for A Vel does not work well for a Ori 
because the terminal velocity of the wind (~10 km s _1 ) lies within the chromo- 
spheric line profile, which is broadened by turbulence of ~20 km s _1 . However, 
another approach is available. Because a Ori has a large angular diameter (the 
largest of any RSG), its wind can be resolved by the Very Large Array (VLA) 
radio telescope. The resolved radio continuum is a useful diagnostic because the 
source function is thermal (S u = B u ) and controlled by local con ditions, which 



greatl y simplifies modelling of the outer atmosphere and wind. lHarper et al 



( 200ll j used multi- wavelength VLA radio continuum fluxes, constrained by ad- 
ditional published r adio and mid-IR flux es, to derive an empirical model of the 
wind of a Ori. The Harper et al. ( 200ll ) density model, assuming a spherical, 
steady flow, gives a mass-loss rate of M = 3 x 1CU 6 M yr _1 and a slow wind 
acceleration with a terminal velocity of ~10 km s _1 . 



3.3. Case Studies of Two Supergiant Binaries: £ Aur & VV Cep 

Observations of RSGs in binary systems provide a unique opportunity to derive 
spatially resolved information about their outer atmospheres. RSG binaries are 
coeval systems with stars of roughly similar mass. The more massive star has 
evolved off the main sequence into the RSG region of the HRD, while the less 
massive component remains on the main sequence as a (typically) B star. In 
the ultraviolet (UV), essentially all the flux comes from the hot companion. 
Since the companion is very much smaller (in physical size, not mass) than 
the RSG, the UV line of sight to the companion samples only a small volume 
of the supergiant's wind. Along sightlines close to the supergiant's limb, a 
"chromospheric" absorption spectrum is seen superimposed upon the typically 
featureless B-type continuum. For the special case of eclipsing binaries, the 
sightlines during ingress and egress probe deep into the supergiant's atmosphere 
and provide an unparalleled look at atmospheric structure on very small spatial 
scales (about 1-2% of the supergiant's radius). 

The potential of these RSG binaries (the £ Aur stars with G-K supergiant 
primaries, and the VV Cep stars with M supergiant primaries) has been rec- 
ognized for decades but never fully realized, largel y because these observations 
were done in the blue- violet region of the spectrum ([Wilson &: Abt|[l954l ; IWrightl 
197dl ). In this spectral region, both stars typically have comparable fluxes, and 
the composite spectrum must be disentangled before useful analysis can be done. 
Also, nearly all of the strong resonance and low-lying atomic transitions, which 
are most useful for mapping the chromosphere, lie in the UV spectrum inacces- 
sible from the ground. The situation has changed dramatically in recent decades 
starting with the launch of IUE in 1978, and continuing with the HST/ GHRS 
and STIS spectrographs, and the recent Far Ultraviolet Spectroscopic Explorer 
(FUSE) mission. These space missions revolutionized stellar astronomy by mak- 
ing the UV spectra of stars accessible. 



3.4. £ Aurigae 

In the mid-1990s, £ Aur (orbital period P = 2.66 yr) was observed at a series of 
11 epochs over 2 orbital periods (Figure [2]) by the HST/GHRS. These observa- 
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Figure 2. Left: The relative orbit of £ Aurigae, as viewed in the orbit plane 
(top), and as projected on the sky (bottom). The line of sight to Earth is 
indicated by the dotted line. Positions of the B-type companion at HST 
observation epochs are shown. Right: The spectrum of £ Aur near 1700 A at 
various orbital phases, as observed by the HST/GHRS echelle spectrograph. 
Epochs 4 and 5 show the chromospheric spectrum seen near eclipse. 



tions permitted the circumstellar density profile to be mapped, the mean wind 
law and velocity structure to be derived (assuming steady, spherical flow), and 
provided definitive evidence of substantial orbit-to-orbit variability in the wind 
structure. This analysis gave a mass loss rate M = 5 x 10 -9 M yr -1 with a 
terminal velocity of 70 km s . The development of the chromospheric spectrum 
as eclipse ingress progresses is also shown in Figure [2j It is evident that line 
profiles show complex structure, especially near eclipse where two components 
are often present in the line core (see Epoch 4 in particular). 



3.5. VV Cephei 

The UV spectrum of VV Cep was observed in total eclipse, for the first time, 
in November 1997 with the newly commissioned HST/STIS echelle spectro- 
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graph. Observations of VV Cep continued on for a total of 21 epochs, end- 
ing just past first quadrature in 2003. In eclipse, VV Cep has a rich emis- 
sion line spectrum with more t han 2000 features between 1300 and 3150 A 
( Bauer. Bennett. Sz Brownl 120071 ). Immediately after egress, VV Cep shows 
a complex chromospheric absorption spectrum, while far from eclipse (near 
quadrature), an inverse P-Cygni "shell" spectrum becomes prominent. In partic- 
ular, the chromospheric lines sometimes show doubling of the line profiles, indi- 
cating the presence of complex velocity structure in the lower chromosphere. De- 
spite this complexity, reliable line-of-sight (tangential) column densities through 
the chromosphere can be obtained from analysis of the prominent damping wings 
of the stronger atomic lines, e.g., Mgn h & k, Hi Lyman-a, and the stronger 
lines of Fe II. These tangential column densities can then be used to construct 
a density model of the chromosphere and wind of VV Cep. The resulting wind 
model has a mass-loss rate of M = 1.5 x 10~ 6 M Q yr _1 with a terminal velocity 
of Voo = 20 km s _1 . 



4. The Wind Velocity Law 

We can use these RSG analyses to piece together an idea of the actual wind 
velocity law. The binary technique excels at observations in the chromosphere 
and the base of the wind, where pertubations due to the companion are minimal. 
But, given the separation of stars is only about a/-R*~ 6 for both ( Aur and 
VV Cep (where a is the semi-major axis, i?* the stellar radius), it is unlikely 
that the binary approach provides much useful information at radial distances 
beyond r ~ a/2 ~ 3i?*. In this regard the radio continuum flux method is 
complementary since (at present) the spatial resolution is limited to about 3i?* 
for a Ori. Similarly, the line profile analysis used for A Vel works best for the 
strong lines with well-defined absorption troughs with well-defined blue edges 
and so is well-suited for recovering structure in the outer wind. 

One can get a sense of the complete wind laws for K and M supergiants 
by combining the results of these analyses in a single diagram (Figure [3]) . The 
dashed curves derived for the binary supergiants are only useful out to a distance 
of about 3-R*, while the solid curves of the single stars are only meaningful 
beyond 2-3 i?*. Presumably combining the useful parts of each of these curves 
should approximate the complete velocity law of K and M supergiants. 



5. Line Doubling in the Chromospheric Spectrum 

Despite these results, we still know little about what drives these supergiant 
winds. Some of the assumptions adopted are manifestly incorrect, e.g., that 
of a steady wind, given the obvi ous line profile chang es in the UV spectrum 
of A Vel between 1990 and 1994 ((Carpenter et all 11993 . Fig. 10). The cores of 



chromospheric a bsorption lines a r e typi cally double (e.g., see the Epoch 4 spec- 
trum, Figure E]) . iMcKellar et al. ( 19591 ) presented the first detailed description 



of RSG chromospheric line doubling, as observed in the Call K-lme of the RSG 
31 Cygni. Line doubling is seen in our HST observations of VV Cep (Figure 0]). 
The profile shape and position vary somewhat from epoch to epoch, but the 
separation of the components remains approximately constant. For VV Cep, 
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Figure 3. A comparison of wind velocity laws for four RSGs. The results for 
single stars (solid) and binaries (dashed) are shown together. K supergiants 
are shown in blue and M supergiants in red. v oo is the terminal velocity. 

averaging over several lines and epochs gives a mean separation of the compo- 
nents of Av = 14.5 ± 2.2 km s . These components straddle the supergiant's 
RV, with the mean blue component at vuue = —11-3 ± 1.8 km s -1 (systemic). 

The relative strength of these components is strongly correlated with lower- 
level excitation potential. Low excitation lines have a dominant blue component, 
but the strength of the red component increases with excitation potential. Lines 
from levels above 1 eV show predominantly red components. The line-doubling 
phenomenon appears to be ubiquitous in RSG chromospheres and needs to be 
understood before more realistic models can be constructed. 
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Figure 4. The variation of chromospheric line profiles in VV Cephei with 
lower-level excitation energy E\ ow . The spectra are labelled by HST epoch 
number on the right. Epoch 10 on 1998 Dec 1 is the first observation fol- 
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each epoch. Wavelengths of atomic transitions are in vacuum. 
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